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Abstract 

We report on the long-term evolution of the spin period of the symbiotic X-ray pulsar 
GX 1+4 and a possible interpretation within a model of quasi-spherical accretion. New 
period measurements from BeppoSAX, INTEGRAL and Fermi observations have been com- 
bined with previously published data from four decades of observations. During the 1970's 
GX 1+4 was spinning up with the fastest rate among the known X-ray pulsars at the time. 
In the mid 1980's it underwent a change during a period of low X-ray flux and started to 
spin down with a rate similar in magnitude to the previous spin up rate. The spin period 
has changed from ~110 s to ^160 s within the last three decades. Our results demonstrate 
that the overall spin down trend continues and is stronger than ever. We compare the ob- 
servations with predictions from a model assuming quasi-spherical accretion from the slow 
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wind of the M giant companion. 



1 Introduction 

Accreting X-ray pulsars are highly magnetized neutron stars in a binary system, accreting 
matter from their companion star. The mass transfer can take place via Roche-Lobe overflow, 
strong stellar winds for giant stars or the Be emission mechanism. These accreting pulsars 
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radiate predominantly in the X-ray band, and the radiation is modulated by the stellar 
rotation of the pulsar. For a review on this subject, see e.g., Nagase [13] . 

GX 1+4 is an accreting X-ray pulsar which was discovered in 1970 by a balloon X- 
ray observation at energies above 15 keV showing pulsations with a period of about two 
minutes (Lewin et al. [12]). It was one of the brightest X-ray sources in the Galactic center 
region. The composite emission spectrum of GX 1+4 indicated that the object was almost 
certainly a binary system, consisting of a symbiotic red giant and a much hotter source 
(Davidsen et al. [3]). Chakrabarty & Roche [2] confirmed the optical companion to be the M 
giant V2116 Oph. 

More recently, observations have shown a well- determined, 1161 day period, single- 
line spectroscopic binary orbit (Hinkle et al. [II]). The inclination of the orbit is, however, 
unknown. Therefore, on the basis of similar characteristics seen in other X-ray binary systems 
the mass of the neutron star is inferred to be ~ 1.35 M@. The mass function of the system, 
combined with the abundances of the M giant star and the assumed mass for the neutron 
star, indicate a mass of ~ 1.2 Mq for the M giant star. This implies that the M giant star 
is a first ascent giant which does not fill its Roche Lobe. In consequence, the neutron star 
is capturing the slow stellar wind of its companion, which makes V2116 Oph quite different 
from other low-mass X-ray binaries (Hinkle et al. [H]). Therein, GX 1+4 is the first low 
mass X-ray binary with an M giant optical companion, and the first symbiotic X-ray pulsar 
ever observed. 



2 Results 

The main result obtained during this work is the plot shown in Fig. [l] which currently is the 
most updated and complete series of pulse periods of GX 1+4. The available measurements 
of the pulse period of the neutron star of GX 1+4 span a time period of about 40 years giving 
us a unique insight in the pulse period evolution of this symbiotic X-ray pulsar. To obtain 
this result we have combined all measurements that we are aware of from the literature with 
the pulse periods determined from BeppoSAX/WFC, INTEGRAL /ISGRI and Fermi/GBM 
data. 

The system has a peculiar spin history. During the 1970's GX 1+4 was spinning up 
with the fastest rate (P ~ 7.23 x 10 -8 s/s) among the known X-ray pulsars at the time (see 
e.g., Nagase [13]). During that period the source had an X-ray luminosity between 9 x 10 38 
and 2.25 x 10 38 erg s" 1 (White et al. p2]) (the luminosity value have been corrected to 
the distance range (3-6 kpc) (Chakrabarty &; Roche [2]) ). However, EXOSAT observations 
in 1983 revealed an extended low state with an X-ray luminosity of < 4 x 10 36 erg s _1 
(Hall & Davelar [TO]). Around the same time GX 1+4 started to spin down with a rate 
similar in magnitude to the previous spin-up rate (see Fig. [T]) and continues spinning down 
up to date. Indeed, the spin down rate measured by Fermi is the strongest spin down rate 
observed to date (P ~ 1.08 x 10~ 7 s/s), and the pulse period of the pulsar has increased by 
about 50% during the last ~ 30 years to the largest value ever known for this source (~ 160 s) 
(Gonzalez-Galan et al. [9]). 
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Figure 1: Long-term pulse period of GX 1+4 (Gonzalez-Galan et al. |9] and references 
therein). Note that the period increases from top to bottom. 
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Figure 2: GX 1+4 data. Top left plot: Top: CGRO/BATSE daily lightcurve in the energy range 20-40 
keV. Middle: pulse periods derived from CGRO/BATSE data. Bottom: Residuals of the periods from a linear fit. Top 
right plot: Top: CGRO/BATSE daily lightcurve in the energy range 20-40 keV. Middle: pulse periods derived from 
BeppoSAX/WFC data. Bottom: Residuals of the periods from a linear fit. Bottom left panel: INTEGRAL /ISGRI 
data. Top: Average flux per pointing in the energy range 18 to 40 keV. Middle: 20-40 keV pulse period. Bottom: 
Residuals from the pulse periods using a linear fit. Bottom right plot: Top: Swift/BAT daily averaged light curve in 
the energy range 15-50 keV. Middle: Pulse periods derived from Fermi/GBM data. Bottom: Residuals from the pulse 
periods using a linear fit. Note that in all these plots the period increases from top to bottom. 
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Apart from the long-term plot, looking at the behavior along shorter time scales (Fig. [2]), 
it is possible to observe that even when the X-ray flux has variations, pulse period evolution 
seems to be almost linear. 

The unusual long-term spin behavior of GX 1+4 has attracted considerable interest for 
many years. Studying the pulse period evolution in an accreting X-ray pulsar and relating 
it to, e.g., the luminosity changes, allows to test models of accretion and to gain insights on 
the physical processes taking place in this kind of systems. In Gonzalez-Galan et al. |9] we 
extend the investigation of the spin period history of GX 1+4 with new observations obtained 
by BeppoSAX, INTEGRAL and Fermi, we correlate X-ray flux variations with pulse period 
variations of the system using the same observations as well as those obtained with Swift and 
BATSE, and provide a likely explanation for the spin history seen. 

3 Accretion models 

Standard accretion theory by Ghosh & Lamb [8] assumes the formation of a prograde disk 
around the neutron star. This prograde disk accelerates the neutron star through the angular 
momentum of the accreted matter. Therefore, spin down of the pulsar is only possible in the 
quasi-equilibrium state, and there is a positive correlation between accretion rate (M) and 
the acceleration of the neutron star. Assuming M oc Fx, the positive correlation i> oc Fx 
should be observed in the data even during spin down. 

Transient disks with an alternating sense of rotation are known to form in numerical 
simulations in binary systems fed from stellar wind (Fryxell & Taam [7]), therefore it is also 
possible the formation of a transient retrograde disk in GX 1+4 as it is a wind-fed source. 
Indeed, a retrograde accretion disk was proposed for the first time by Makishima et al. [13] 
to explain the spin down observed in the neutron star of GX 1+4. This model assumes a ret- 
rograde disk around the neutron star decelerating the pulsar through the angular momentum 
of accreted matter, therein, a negative correlation — z> oc Fx should be observed in the data. 

There is a third model proposed for the first time to explain the behavior of GX 1+4 in 
Gonzalez-Galan et al. [9], the quasi-spherical accretion model. This model predicts a quasi- 
static atmosphere around the neutron star (instead of a disk) which carries away angular 
momentum from the neutron star magnetosphere but allows accretion (in contrast to the so- 
called subsonic propeller regime (Davies & Pringle[l])); and a negative correlation (— v oc Fx) 
during spin down (see Shakura et al. [T7]; Postnov et al. |16j). The formation of this kind of 
"atmosphere" is only possible for wind- fed sources like GX 1+4. 

4 Discussion 

The correlations between X-ray flux and pulse period evolution predicted by the models 
proposed in Section [3] can be tested by comparison with a regular sequence of X-ray fluxes 
and pulse periods which can be found in BATSE and Fermi data (see Fig. [2]) for GX 1+4. 

A cross-correlation of v obtained from Fermi/GBM data and the contemporaneous 
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Figure 3: The discrete correlation function between u obtained from Fermi/ GBM GX 1+4 
data and the Swift/BAT 15 - 50 keV X-ray flux shown in bottom right panel of Fig. [2j The 
minimum near zero lag implies a negative correlation between the pulse frequency derivative 
and the X-ray flux (Gonzalez-Galan et al. [9]). 

X-ray flux measured by Swift/BAT finds a minimum near zero lag (see Fig. [3]) implying a 
negative correlation between this two parameters for GX 1+4. Furthermore, this negative 
correlation has already been found previously in BATSE data (i.e. Chakrabarty et al. [1]). 
This negative correlation discards the standard accretion disk model as a possible interpre- 
tation of GX 1+4 behavior. 

Due to several reasons, including this negative correlation, different authors have tried 
to explain GX 1+4 behavior assuming the formation of a retrograde disk around the neu- 
tron star by matter captured from stellar wind of the secondary M-giant star in GX 1+4 
(e.g. Makishima et al. [13], Dotani et al. (5], Chakrabarty et al. pQ, Nelson et al. [15], Fer- 
rigno et al. [6], etc.). The weak points of the retrograde disk interpretation, however, include 
the puzzling long-term (~ 30 years) stability of such a retrograde disk. 

This puzzling long-term stability of the retrograde disk interpretation has led us to 
a new model, the quasi-spherical accretion model, which is able to explain the negative 
correlation found in the data without the formation a retrograde accretion disk (Gonzalez- 
Galan et al. [9]). 

5 Conclusions 

New measurements by INTEGRAL and Fermi confirm the overall spin down of the neutron 
star in GX 1+4 observed over the last decades. This spin down is stronger than ever ob- 
served before with an increasement of the spin period of ~ 50% during the last ~ 30 years 
(Gonzalez-Galan et al. [9]). The analysis of Fermi data reveals a negative correlation be- 
tween spin frequency derivative and the X-ray flux simultaneously measured by Swift/BAT 
monitor (Gonzalez-Galan et al. [9]), which confirms the correlation of instantaneous spin- 
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down torque with X-ray flux discovered by BATSE (Chakrabarty et al. pQ) discarding clearly 
the standard accretion disk model for this system. We have shown that in GX 1+4 not 
only the retrograde disk accretion is a possible explanation but the quasi-spherical accretion 
onto the neutron star from the stellar wind of the M giant companion is likely to take place 
(Gonzalez-Galan et al. [9]). 
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